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We report on recent study [T] of the role of the kinematical constraint 
in the CCFM equation and its non-linear extension. We compare numerical 
results obtained by solving the CCFM equation and argue that kinematical 
constraint represents an important correction. 

PACS numbers: PACS numbers come here 

1. Introduction 

In the regime of high energy scattering, where the center of mass energy 
is larger than any other available scale, perturbative approach to processes 
with high momentum transfer allows factorization of the cross section into a 
hard matrix element with initial off-shell gluons and an unintegrated gluon 
density mm- The unintegrated gluon density is a function of the longitu¬ 
dinal momentum fraction x and transverse momentum kx of a gluon. After 
taking into account formally subleading corrections coming from coherence 
of gluon emissions one is lead to the CCFM equation which introduces gluon 
density dependent on hard scale related to probe. In principle, it is an equa¬ 
tion that should be the ideal framework for application to final states at high 
energies and covering DGLAP and BFKL domain in gluon channel. It has 
been implemented in the Monte Carlo event generator [2]. However, so far 
good agreement with high precision data has been successfully achieved only 
in rather inclusive processes like F 2 and Drell-Yan [5]. It is known that on 
the theory side the CCFM physics is still to be completed. For instance: 
the impact of the kinematical effects introducing energy conservation in 
the CCFM evolution have been not investigated in all detail. As it turns 
out from our study it is necessary to revisit the inclusion of the so called 
kinematical constraint [5] into CCFM equation. Only recently the CCFM 
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has been promoted to non-linear equation CHS El therefore allowing for the 
possibility to investigate interplay of coherence effects and saturation [2]. In 
particular the important question is what is the role of the angular ordering 
and kinematical effect in the evolution at the non-linear level. The optimal 
form of initial conditions is also not known. 

We present here a continuation of the work done in HD El. 


1.1. The CCFM equation and the kinematical constraint 
The CCFM equation for gluon density reads: 

l_«o 

fd 2 q f ? 

A(x,k T ,p) = A 0 (x,k T ,p) + as J — J dz 0(p - zq) V(z,kT,q)* 

X 

A s (p,zq) A{x/z,k T ',q), 

where k 1 = |k + (1 — z)q| and the modulus of two dimensional vectors 
transversal to the collision plane are denoted |k| = kx, |q| = qx and x 
is gluon’s longitudinal momentum fraction and a s = N c a s /Tr. Also the 
rescaled momentum is introduced as q = |q| = qr/iX ~ z). 

Where T (z,kx,q) is the gluon splitting function. The function A s(p,zq) 
is the Sudakov form-factor. The non-Sudakov form-factor A Ns{z,kT,q) 
regularizes 1/z singularity. 


The kinematical constraint 

The integration over q in the equation (jT]) , although being constrained 
from below by the soft cut-off Qq, is not constrained by an upper limit 
thus violating the energy-momentum conservation. Moreover in the low x 
formalism one requires that in the denominator of the off-shell gluon prop¬ 
agator one keeps terms that obey |£r’| = kj,. In order to be consistent 
the non-Sudakov form-factor should be accompanied by a kinematical con¬ 
strained limiting the above integration over q. In approximated form it reads 
k^ > zq 2 and at z <C 1 guaranties that \k 2 \ ~ k^- In El ^ has been ex¬ 
tended to region including also the case when z ~ 1. Careful derivation leads 
to following form of the kinematical constraint kj, > * 1 g 5, = z(l — z) q 2 . 
The lower bound on z > x results in the upper bound on qj, < k^/x ~ s 
providing local condition for energy-momentum conservation. We include 
the kinematical constraint in the CCFM equation. The non-Sudakov form- 
factor after inclusion of the full form of the kinematical constraint assumes 
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form 

1 (&_ n ( (1 ~ z')k\ _ , 
z' \(l-z) 2 q 2 
r <i a ' 2 1 W 

j ^e(k 2 T - q ’ 2 )@(q'-z'q)\ . 

Please note the presence of the function Z )‘ The authors of [Bj 

solve the CCFM at small z limit therefore the function 0 ^ ~ z ^j 

neglected and in most of the phenomenological and theoretical applications 
of the CCFM this term is neglected 0 [TOl H3 [HI [T5j [161 El EE]- The 
following form of non-Sudakov form factor is usually being used: 

A Ns(z,k T ,q) = exp ^-a s ^ j ^-©(fc 2 - q' 2 )@{q' - z'q^j (3) 

= exp a s log (^) log 

where zo = 1 ; if (kr/q) > 1; Zq = kr/q, if z < {kx/q) < 1; Zq = z , if 
( kr/q ) = z. The discussed above 0-function is not taken into account. 



Ans(z, kr,q) = exp \ - a s 



Saturation effects and kinematical constraint combined 

To account for gluon recombination at large gluon densities the CCFM 
equation has been promoted to non-linear equation by including a quadratic 
term BUIE! and it reads: 


Q o 


(1 — z)q 2 


— - z )e(p - zq) 


A(x,hr,p) = A 0 (x,k T ,p) +a s J 
V(z,k T ,q) As(p, zq) (a (^,k T ',qj ~ 5 (q 2 - }%) q 2 A 2 , (4) 


d z q 
7rq 2 


where kp = kp/{ 1 — z) and we included the kinematical constraint of the 
form © in the kernel. Simpler versions of the equation above have been 
already analyzed in pj2| and it has been observed, that the equation leads 
to phenomenon called saturation at the saturation scale nans] and the 
saturation strongly suppresses the gluon density at low x and low kp- The 
natural question arises how are these results modified when some of the 
approximations are not taken and how are they modified if the kinematical 
effect is imposed in the full form. 
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Fig. 1. Relative ratio of CCFM and KGBJS solutions. Distributions with definite 
p for varying value of kp. 
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Fig. 2. Relative ratio of CCFM and KGBJS solutions. Distributions with definite 
p for varying value of kp. 


2. Numerical results 

Linear equations 

We use an initial condition which includes resummed virtual and unre¬ 
solved contributions, according to [12] and [1TJ, in the form Ao ( x,kr,P ) = 
A Ar(x, hp) A s(p, Qo)/kp, with A = 1/2 and Ar(j z, kp) is the Regge form- 
factor. 

The observation we make from the plots (like plots in Fig. [1|) is that the 
solutions of equations we study differ significantly. The solutions exhibit also 
similar features. Solutions of both versions of the kernel with kinematical 
constraint exhibit a local maximum as functions of kp and p. The positions 
of local maxima in the plots of p dependence are correlated with the value of 
kp, with a shift to higher kp for the solution of the equation with kinematical 
constrain ^-function included. The peak can be explained by the fact that 
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the contribution of the integral on the right hand side of the equaiton peaks 
at around kp ~ p■ The peak therefore is a result of presence of 6 (p — zq ) 
- angular ordering condition. Similar peaks are present also in the plots of 
kp dependence and resemble Sudakov suppression of kp scales of the order 
of p m • However, in the case of without the ^-function it seems, that the 
position of the peak does not depend on the value of p. The peak observed in 
solution of the equation without kinematical constrain 0-function is ’hidden’ 
under the result of the evolution. We can conclude that the peak in the kp 
dependence is an interplay result of inclusion of the explicit kinematical 
contraint via 0-function factor and the Sudakov effect. 


Non-linear equation 

We set the parameter characterizing the strength of the non-linear term 
R the value R = \/1/f GeV in the equation Q. By comparing the CCFM 
and KGBJS equations we see, that the kinematical constraint suppresses 
the growth of the gluon so much that the non-linear effects enter only at 
very low x. Observations made in previous paragraphs are confirmed in 
2-dimensional plots (Fig. [2]), where we plot absolute relative difference of 
two amplitudes, solutions of the CCFM and KGBJS equations, defined by 
the quantity 


/3 (x, k T ,p) 


I ■Accfm ( x , kp,p) - Akgbjs {x, 
■A-ccfm ( x , kr,p) 


( 5 ) 


The function (5 ( x , kp,p), introduced before in [T2j, can be used to measure 
the strength of the non-linear effects and to define a saturation scale using 
the conditions: 


)3 (x,Q s (x,p) ,p) = const., f3 (x,kp, P s ) = const. (6) 

The second condition in Q defines p-related saturation scale. 

The conditions above can be seen as equipotential lines in the plots in 
Fig. El where different equipotential lines correspond to different constants 
on the right-hand side of the equation above. The change in the slope of 
the (3 (x, Q s ( x,p ) ,p) at around kp = p reported in jTS], can be seen clearly 
in Fig. El can be understood in the context of the peak at p ~ kp (Fig. (U). 

By comparing the plots Fig. El to analogous plots in [12] we see that 
the main features are very similar. We therefore conclude that the low — x 
approximation of the KGBJS and CCFM equations taken in [T2] does not, at 
least, modify the relative difference between linear and non-linear equation. 
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